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We review the basic properties of the Haldane phase in spin-1 Heisenberg anti ferro¬ 
magnetic chains, including its persistence in quasi-one-dimensional geometries. Using 
large-scale numerical simulations, we map out the phase* diagram for a realistic model 
applicable to experimental Haklane compounds. We also investigate tlie effect of dif¬ 
ferent chain coupling geometries and confirm a general mean field universality of the 
critical coupling times the coordination number of the lattice. Inspired by the recent 
development of character ization of symmetry protected topological states, of which the 
Haldane phase of spin-1 Heisenberg antiferromagnetic chain is a preeminent example, 
we provide direct evidence that the quasi-one-dimensional Haldane phase is indeed a 
non-trivial symmetry protected topological state. 

Keywords: Low dimensional magnetism; Haldane conjecture: quantum Monte Carlo. 

1. Introduction 

Low dimensional interacting spin systems have long been a laboratory for the dis¬ 
covery of novel quantum states of matter. Enhanced quantum fluctuations due to 
reduced dimensionality enable the appearance of multiple quantum phases with 
unique characteristics driven by the interplay between strong interactions, ex¬ 
ternal (e.g.. pressure and magnetic field) and internal (c.g., crystal field effects) 
potentials and lattice geometry that are suppressed in higher dimensions. The 
relative simplicity of the microscopic models facilitates the development of exact 
analytical solutions in many cases and powerful field theoretic and computational 
approaches in others providing greater insight into the emergence of these complex 
phases and their physical properties. Concurrent rapid advances in the synthesis 
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mid characterization of low-dimensional (|inmtuin magnets have kept the study of 
these systems one of the most active frontiers in condensed matter physics (see 
Landee and TnrnbulP for a recent pedagogical review). One of the most remark¬ 
able results in the study of quantum spin models one that substantially enhanced 
our understanding of long range order in quantum many body systems is the 
pioneering work by Haldane03 By studying the non-linear sigma model in (1 + 1) 
dimensions. Haldane conjectured that the ground state of the one-dimensional (ID) 
Heisenberg ant iferromagnet (HAFM) has gapless excitations for half-odd integer 
spins, whereas that for integer spins Is separated from all excited states by a finite 
spin gap (Haldane gap). Haldane's conjecture lias inspired numerous theoretical 
studies of integer spins in low dimensions, primarily S = 1 spins where the Haldane 
phase is most robust. These include chain mean field theory (CMFT) SB exact 
diagonalization!3H density matrix renormalization group (DMRC!){^ 0 - ‘ -1-113 and 
quantum Monte Carlo (Q.MC) simulationsH0E3®lII3 

The theoretical studies have been complemented by the discovery of 
several quasi-one-dimensional (QlD) spin S = 1 quantum magnets, 

such as AgVP 2 S„l M 2ll F ,| F 2 l NDMAl'P NENB,^ N'ENP0 ,, P 1-‘ F3 N1NOP 
PbNi 2 VaOsP?3 SrNiiViOsF™^ 1 TMNINand Y 2 BaNi0 5 M In these ma¬ 
terials, the magnetic ions are arranged in chains, with weak but finite inler-chain 
couplings, which affect the ground state phases. Additioually, in most known S — 1 
magnets, the ubiquitous Heisenberg exchange is complemented by a single-ion 
anisotropy. The expanded Hilbert space of the S = 1 spins, and the interplay 
between multiple competing interactions, external magnetic field and different lat¬ 
tice geometries result in a rich variety of ground state phases. In addition to the 
gapped Haldane phase, examples of exotic quant um states realized in low dimen¬ 
sional interacting spin systems include experimentally realized Bose- Einstein Con¬ 
densation (BEC) of magnonsP 6 !** 7 ! quantum par.uuagnotpH : * <> i 4 lll and the recently 
proposed spin supersolidHIB and fcrroucniatitST phases. Recent advances in syn¬ 
thesis techniques have made it possible to engineer quasi-low-dimcnsioual materials 
where the “effective dimensionality" (that is, inter-chain or inter-layer couplings) 
and Hamiltonian parameters (such as the ratio of exchange interaction and single¬ 
ion anisotropy) can be coiitroHedl J4 l 4 'l This raises the possibility of preparing ma¬ 
terials with desired predetermined properties. The search for such tailor-made ma¬ 
terials has grown in recent, years, as these are believed to drive the next generation 
of electronics. In addition to condensed matter systems, rapid advances in the field 
of ultracold atoms in optical lattices have opened up a new front ier in the study of 
interacting many-body systems in arbitrary dimensions. The unprecedented control 
over number of atoms, interactions, and lattice geometry mokes it iui ideal testbed 
for preparing and studying novel quantum states. 

The most remarkable property of the Haldane phase is that even though there 
is no long-range magnetic order all spiu-spin correlations decay exponentially 
there exists a hidden or nonlocal order measured by the string order parameter 
introduced by den Nijs and Roiumelst®9. The introduction of a nonlocal order 
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Fig. 1. Illustration ol the AKLT state for a six site open chain. Large .shaileil circle represent 
spin 5 = 1 degrees of freedom ccuu|Mised of the symmetric combination of two spin 5=1/2 spins 
(small solid circles). In the AKLT state, pairs of S = 1/2 spins on neighboring sites form a singlet 
state (blue lines). For open boundaries, this leaves an unpaired 5 = 1/2 degree of freedom at each 
boundary. In the thermodynamic limit, these spins are completely decoupled from the rest of the 
lattice, leading to a fourfold degenerate ground state. 


parameter was a marked departure from the conventional practice of characterizing 
quantum many body states in terms of local order parameters and ushered in the 
study of topological phases an area of intense current research. Much insight into 
the nature of the Haldane phase can be gained from the AKLT (AiHeck-Kennedy- 
Licb-Thsaki) state the exact ground state of a ID S' = 1 spin chain where the 
near-neighbor Heisenberg interaction is supplemented by an additional interaction 
between neighboring spins.23 The AKLT state may be understood by noting that 
the spin at each site can be thought of as a symmetric combination of two .9 = 
1/2 spins. Pairs of these spins on neighboring sites form a singlet on each bond. 
For a periodic chain, this forms a unique valence bond solid ground state a 
singlet, on each bond with a gap to lowest excitations. Hut for open chains, there 
remains ,ui unpaired .9 =1/2 moment at each end which Is doubly degenerate (see 
Fig. [TJfor an illustration of this state). Consequently, the gromid state is gapped in 
the “bulk" iuul lias degenerate gapless edge states. The state Is protected by this 
Z 2 x Z 2 symmetry. Using the current terminology of topological states, this state Is 
asymmetry protected topological (SPT) phase. The Haldane phase is adlabaticnlly 
connected to t he AKLT st ate in other words, it has the same qualitative character. 
In fact, the Haldane phase is widely recognized as the earliest and best understood 
SPT phase in interacting many hotly systemsli£li2] 

What happens to the Haldane phase in the presence of additional interactions, 
such as. inter-chain coupling and single-ion anisotropy? The question is not simply 
of academic interest all real materials have nonzero interchain coupling and in 
many quantum magnets, the crystal elect ric field lifts the degeneracy of the lo¬ 
cal Hilbert space in the form of a single-ion anisotropy. Theoretical studies have 
shown that both these interactions destroy the Haldane phase at sufficiently strong 
interaction strengths the single-ion anisotropy drives a transition to a quantum 
paramagnetic phase whereas interchain coupling favors long-range AFM ordering 
but the ground state remains gapped up to finite values of the couplings. What is 
the nature of the gapped ground state away from the isotropic Heisenberg point? 
Does it still retain its SPT character? The string order parameter the only defini¬ 
tive probe for the Haldane phase is strictly defined only in one dimension and 
its extension to coupled chains cannot be trusted without corroboration from addi- 
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lional measurements. Fortunately, recent advancements in the study of SPT phases 
have yielded new probes that provide Letter insight into the nature of the putative 
Haldane phase away from the isotropic Heisenberg point. In tills brief review, we 
summarize some of our recent progress towards fully characterizing the Haldane 
phase in quosi-one-dimensional geometries. 

2. Model 

We consider a model of low dimensional quantum magnets consisting of an array 
of weakly coupled spin-1 HAFM chains. Tliis model can be described by a spa¬ 
tially anisotropic Hamiltonian of S = 1 spin operators with nearest-neighbor spin 
exchange interactions and on-site single-ion anisotropy: 

n = j'£Z-s J + j''£Si-§i + DY,(S‘)\ (i) 

<ij> 101 ' 

The first two sums are over all pairs of interacting spins ij, with ( ij) and (ij) refer¬ 
ring to nearest neighbor spin pairs along and perpendicular to the chain direction, 
respectively, while the third sum is over all lattice sites i. Here, we limit ourselves 
to nearest neighbor interactions on bipartite lattices to avoid the sign problem in 
QMC, and tune the interactions ./ and ./’ to generate systems of weakly interacting 
spin chains. This is illustrated in Fig. [2] for a spatially anisotropic square lattice 
composed of weakly interacting spin chains. In this paper, we consider several bipar¬ 
tite lattices in the quiisi-one-dimensional limit. For each lattice, we define a chain 
length L as the system size along tin* chain direction, and a perpendicular length 
Lx. such that for 2D lattices the total number of spins N = L x L± while for 3D 
lattices N = L x L± x L±. The ratio of these two length scales forms the aspect 
ratio It = L/L±, which we set at ft = -I unless otherwise stated. We henceforth set 
the intraclmin iuteractlou strength to unity (J = 1), and study the ground state 
properties of the above model for various single-ion anisotropy D and interchain 
coupling 

3. Methods 

Quantum Monte Carlo methods map quantum systems onto equivalent classical 
representations in d + 1 dimensions, where the configuration space can be stochas¬ 
tically sampled by Markov chain Monte Carlo methods.^ The numerical results 
presented here were obtained using the stochastic series expansion (SSE) formula¬ 
tion of Q.MCl 51 ! 52 ! Briefly, the density matrix can be expanded as a Taylor series, 
t~ (— ffH)' /l\. For a spin Hamiltonian consisting of onsite fields ami |mir 

interactions, we can write -H — where the stun is over all pair interac¬ 

tions. or bonds. The powers of V appearing in the Taylor expansion of c -JM can 
then be written as a sum over products of bond operators, (~7t )' = J~[ ( Hb{i,c)- 

where C represents a sequence of bond operators called the operator shiny. Note 
that t he overall sign of a particular operator string depends on the number of bond 
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Fift. 2. Illiberal luii of a spatially anisotropic square lattice composed of spin chains with dominant 
intrachuiu coupling J (l»lu«? lines) uml weak interdiuin cuupliug J' (mil Hues). 


operators with negative weight. To avoid the sign problem in QMC. all the bond 
operators must therefore be positive definite, which in practice limits the study of 
Heisenberg autiferroinagnets to bipartite lattices, where a suitable sublattice ro¬ 
tation is available. The configuration space consists of all possible rearrangements 
of the operator string, the length of which takes into account the expansion vari¬ 
able t. while its sequence dictates t he states a appearing in the partition function. 
Z — y {ole* 1 ™ a). Loop algorithms exist to efficiently sample this configuration 

spaccFWq 

In addition to the SSE method, we implement n projective QMC' method that 
allows us to access both grotmd state expectation values as well as wave function 
overlaps between the ground state and a trivial product state. This method is quite 
similar to SSE. and details are given eLsewhercP^D 


4. Results 

Here we present recent QMC results on the ground state properties of the spin 
S = 1 HAEM described by the Hamiltonian hi Eq. iJTJ. First, we show how finite 
size scaling can be used to accurately determine the ground state phase boundaries 
between the Haldane phase and neighboring magnetically ordered phases. Next, we 
examine the effect uf lattice geometry on the critical coupling of bipartite lattices. 
The phase diagram is then presented in the quasi-one-dimensional regime using a 
sputially anisotropic rectilinear lat tice. This phase diagram is directly compaml to 
the experimentally determined parameters of Haldane gap materials. We also show 
the magnetization curves of the Q1D model for the various zero-field ground state 
phases, mid calculate the low-lying excitation spectrum of the Haldane phase near 
the Haldane gap minimum. Finally, the Haldane phase in QlD geometries is shown 
to belong to the new class of symmetry protected topological states. 
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4.1. Finite Size Sculitiy 

Even with a simulation method such its sign-free quantum Monte Carlo that scales 
polynomials in system size, it is impossible to directly access results for macro¬ 
scopic system sizes. To overcome this difficulty, we employ finite-size scaling (FSS) 
methods to extrapolate our results to the ground state thermodynamic limit. Since 
quantum systems map onto classical systems in d + z dimensions, where : Is the 
dynamic critical exponent, similar procedures can be used to access the ,3 —> x and 
L -» oc limits. For ground states that break a continuous spin symmetry, the spin 
stiffness is a useful observable to determine the boundary of such a gapless phase 
(as per the Goldstoue theorem) with a gapped state of any order. This is because 
the spin stiffness must be finite in the gapless phase and zero in the gapped phase. 
The scaling function of the spin stiffness in the vicinity of a quantum critical point 
Ls known to be 

P.W.S.i) = Z, 2 - J -=/ 1 [p/L\ (J-«.)£■"), (2) 

where g is a Hamiltonian |tarauicter that can tune the ground state system to the 
quantum critical point at y c , v is the critical exponent governing the growth of the 
(spatial) correlation length £ near y c . and : relates £ to the correlation length in 
imaginary time by the relation r =£ : . 

By taking the inverse temper at me to scale as the system size L (with z = 1 for 
the quantum phase transitions considered in this work), it is clear that. ~~ Ls 

independent of system size at the critical point. This Ls demonstrated in the right 
panels of Fig. [jJ where we plot p,L 2 for a 3 + 1 dimensional quantum system and 
use 3 = L to reach the ground state limit. The critical point Ls clearly visible at 
D,- — 0.092(2) where the quautity p%L 2 is independent of system size. 

A similar scaling function is known for the square of the staggered magnetization 
along the ; axis m~ that acts ils an order parameter for the Ling antifcrromagnctic 
phase. It is 

(/}/£■, ( 9 - 9 «)i ,/ -) . (3) 

where r/ Ls the anomalous dimensionality. Thus we can aLso see that 
will dLsplay a crossing behavior for varying L at the phase boundaries of the king 
AFM phase. This Ls shown in the left panels of Fig. [d] for the king AFM to Haldane 
phase transition. 


4.2. Effect of geometry 

Starting from a system of uncoupled chains in the Haldane phase and slowly turning 
on the interchain coupling the system will remain in the gapped Haldane phase 
until a critical coupling ./' is reached at which point the gap closes ami long-range 
magnetic order develops. Within the chain mean field approximation (CMFA), the 
critical coupling .// depends only on the coordination number of chains n. Thus, 
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Fig. 3. Fiuite-iuc Tilling of the Muggi'icil uuiguctbuitkin (mj, bin. pane In) and spin stilTuts- (/»,, 
right panels) a! tin' bang AFM to Haldane uml Hulilaui' tu XY AFM pliase transitions, re¬ 
spectively. Interchain coupling is Axed at J' = 0.012. The critical points D; = — 0.04(1) and 
D' u = 0.092(2) are determined by the crossing criterion (upper panels) assuming a dynamic crit- 
iral exponent r = ! Good finite-sin.' stulmg collapse is achieved near the critical points ussuining 
mean field exponents (»/ = 0 and i» = 1/2) fur a cuutiuuuus phase tiansitkm at the upper critical 
dimension d+ : = 1 (lower panels). Inverse temperature is scaled as 3 = L to ensure convergence 
to the giounil state limit. 


the quantity nJ' might be expected to be universal for unfrustrated lattices. In 
Fig. [j| we show results for chuius arranged into square (n = -I) and honeycomb 
(;t = 3) superlatticcs. While they give different results for the critical coupling 
(./' = 0.0162(4) and ./' = 0.0229(0). respectively), the sealed values nj' c are almost 
in statistical agreement: 0.0648(16) vs. 0.0687(18). These can he compared to the 
CMFA vjiltte nJ' f as 0.051 that nets us a lower hound.® 


4.3. Phu.se Diagram 

The ground state phase diagram of P is shown in Fig. [o] For small |D| and ./' 
the system is in the Haldane phase. For sufficiently strong interchain couplings, the 
Haldane gap is quenched and three-dimensional long-range magnetic order sets in. 
This magnetic order is the Neel antiforroinagnotic stat e in the ease of isotropic spins 
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Fig. 4. Finite-size scaling of tin* spiii stillness (/>.) at tin* Haldane to Neel phuse 1 lansitiun for 
chain coordination (a) n = 1 (i.e. square superlnttice uf dutins). anil (b) n = 3 |i->' hutieycutnb 
superlattice of chains). Single-ion anisolmpy is fixed at the Isotropic point, D = 0. Inverse tem¬ 
perature is scaled as 3 — L/2 to ensure convergence to the ground state limit. Adapted from K. 
Wietsdiem and i*. Seng up t a. JPS Couf. Proe. 3 (2014) 012005. 


(D = 0). while axial (£> < 0) and planar (D > II) anisotropy lead to Lsiug AFM 
and A’l’ AFM states, respectively. Additionally, there is a quantum paramagnetic 
(QPM) phase for large D > 1. The Haldane to XV AFM phase boundary is deter¬ 
mined by FSS of the spin stiffness, while the Haldane to king AFM phase boundary 
Is determined by FSS of the staggered magnetization. In the ease of the Haldane 
to Neel phase transition at the isotropic point (D = 0), FSS of both the spin stiff¬ 
ness and staggered magnetization yield values in agreement up to the statistical 
uncertainty given. 

One feature of interest in the phase diagram of Fig. [a] is the ability to change 
the effective Hamiltonian |taruuictcrs of a given material by the application of hy¬ 
drostatic pressure. This has been demonstrated for NENP by Zaliznyak tt 
However, as can be seen in Fig. [5] NENP under pressure Is actually closer to an 
ideal Heisenberg chain, i.e. it moves towards the origin ami away from the bound¬ 
aries of the Haldane phase. Thus, it is not a good candidate for the observation 
of pressure-induced quantmn criticality, as has been found hi the spin- 1/2 dimer 
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Fin. S. Phase diagram m the D J' (dune with phase boundaries as indicated by the dottiil black 
lines. The borders of the Haldane and QPM ph.-i.-as are obtained as polynomial fits to the present 
work and represent guides lor the eye. Data points for the present work are determined by QMC 
simulations at constant D anil V (solid rod circles and blue squares, respectively). Several Ilnldune 

(iiuipoundi uie plotted us large ciusskntilied symlMils using estimates for D and ./' fium the 
indicated sources. Adapted from K VVieischem and P. Sengupto. Phvs. Rev. Lett. 112. ‘247203 
(2014). 


compound TIC 11 CI 3 J® By contrast, the isostructural coin)Miuiids PbNi.V 2 Os and 
SrNi-jVaOfi are found to already lie very near to a quantum critical point, making 
them prime candidates for the study of pressure-induced quantum criticality. 

In addition to pressure, chemical substitution represents another mechanism 
by wliich the effective Hamiltonian parameters of magnetic compotmds c.ui be 
modified. Tills is already somewhat apparent in the Lsostructural compounds 
PbNi 2 V 2 0a and SrNi 2 V 2 0 8 , where the differing influence of the Pb a+ and Sr 3 * 
ions leads SrNfeVaOs to be closer to the Isiug AFM phase than PbNbV'Os. In¬ 
deed. SrNi 2 V 2 O s was initially thought to magnetically order below T\ — 7/v based 
on experiments on powder samples whereas recent results on poly tryst allim^ 
and single crystal^)have established its low temperature magnetic behavior to lie 
consistent with a non-magnetic Haldane ground state. 

The prospect of chemical fine-tuning Is even more exciting for molecule-based 
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magnets, where the choice of bridging ligands hits been shown to allow for the 
synthesis of a wide rauge of low dimensional structures. Of particular interest Is 
the recently synthesized compound [Ni(HF 2 )(3-Clpy).i|BF.i that is believed to lie 
near the ID Gaussian critical point t'" 1 ' We note here that this critical point is 
well-described by the Tomounga-Luttinger liquid theory with Luttinger parameter 
A' as 1.321 PI 

It is important to take the crystal symmetry into consideration, as it determines 
which symmetry breaking and conserving phases can occur. For example, NENP 
is known to have a rhomlrohedral crystal electric Held component E/J as 0.01 that 
explicitly breaks the onsite f/(l) spin symmetry of H down to Z 2 . Additionally, 
application of a uniform longitudinal magnetic Held Is known to induce a staggered 
transverse field at the Nickel sites due to a zigzag chain structured Also, note that 
the related compound N'ENC falls into the QPM phase since ./ Ddbut because 
E > ./ we do not include it here. 

We also have calculated the phase boundary between the .YV AFM and QPM 
phases from the QlD limit all the way to the spatially isotropic 3D lattice. This 
is shown in Fig. [H] where we also plot the location of some quantum magnets in 
the phase diagram. Note that the phase botmdary Iretween the .YV AFM and 
QPM phases is a line of critical points belonging to the 0(2) universality class in 
3+1 dimensions. Within the QPM phase, application of a magnetic Held leads to a 
field induced quantum phase transition into the canted XY AFM phase. This field 
induced transition has been studied by Zhang et aid and shown to be a mean field 
transition because the effective dimension is V — 3 + 2 (the dynamical exponent is 
2=2 because these field induced transitions belong to the BEC universality (lasi^ZJ). 

The theoretical effect of uniaxial hydrostatic pressure applied to DTN has beeu 
examined, with encouraging signs that DTN may be tuned to a quantum critical 
point (QCP).®3This pressure-induced QCP would be of the 0(2) universality class 
in 3+1 dimensions, as opposed to the field-induced QCP in DTN that has been 
shown to belong to the BEC universality class in 3+2 dimensions.^ There also are 
several compounds much deeper in the QPM phase, such as Sr. |NiPtO ( ,P^ 

The parameters for the compounds shown in Figs, [o] mid [fi] are listed in Ta- 
ble[T] For completeness, we also mention the Haldane chain compounds NDM.AZP^ 
NENBP1 NINAZj® and NINO,® which could also fit into the phase diagram in 
Fig.[5]ucar the related compounds NDMAP .uni NENP but are not. shown here in 
the int erest of clarity. Additionally, the hexagonal compounds CsNiCl, and TMNIN 
are outside the scope of the present model as they cannot, be well approximated by 
a cubic lattice. TMNIN is a Haldane gap material while Cs.NiClj exhibits Haldane 
behavior above its magnetic ordering temperature T,\ — -1.8(1 A'?- 1 Finally, the low 
temperature phase of Tl-Rn-O? has been interpreted as a Haldane chain system 
that spontaneously arises from a fundamentally three-dimensional crystal structure 
due to the orbital ordering.® 
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Fin. G. Phase diagram in the large-/) regime of the D .7' plane. Data puiuts fur the phxise boundary 
are extracted by fiulte-ai&e scaling analysis ul QMC simulations performed at cunsUml J’. Dotted 
line is obtained ns a polynomial lit to tin- QMC data and delineates the XV AFM and QPM 
ph.Lse boundary. 


4.4. Magnetization 

The ground state magnetization curve can he a good way to characterize quan¬ 
tum umgucts tit low temperatures. The direct application of an external magnetic 
field Is tut efficient probe of the underlying magnetic phases of a magnetic com¬ 
pound. For longitudinal magnetic fields, the Hamiltonian in Eq.i[l]i comes with a 
conservation law such that the ground state at any field must occur within a single 
magnetization sector defined by the quantum number M — S;. This can lead 

to several interesting features. First, note that the fully polarized state | | • • • t) is 
.ui exact eigenstate of the Hamiltonian, since the off-diagonal terms come in pairs 
SfSj +S?S% ~ i (.5/5* + .5*5/) such that | f •• • t> isalwaysannihilated,leaving 
only the action of the diagonal terms to set its energy. Additionally, the saturation 
field /?, can often be obtained by setting the energy difference between | T *' * T) and 
the lowest energy state in the M = NS- 1 sector to zero. For the model considered 
here on an anisotropic cubic lattice, it can he shown that B, = 8 ./ + 47' + D as 
long as D is not strongly easy-axis. It has previously been demonstrated that for 
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Table 1. Parameter valuer, fur compound-. shown in Figs. (^] ami [HJ When is nut given in the 
reference, we leave the field ./'/•/ blank and plot the compound as a purely one dimensional chain. 
F«i! NDMAP, we take J' as the average of t he interchain coupling along perpendiculur directions, 
while for PbN^VgOg and SrNi. V 2 0. s vve take as the Ising component of the iuterciutin coupling 
and divide by a factor of 2 lo account for differences in the chain coordination number. Values for 
SraNiPtOc an- averages of results from specific heat and susceptibility data. 


Compound 

Bravais Lattice 

D/J 

J'/J 

J 

Ref. 

AgVP 2 S<i 

monoclinic 

0.000 

1 x 10"* 

780K 

[22] 

NDMAP 

orthorhombic 

0.25 

3.8 x 10" 4 

33. IK 


NENP 

orthorhombic 

0.16 

3 x It)-* 

46K 

[25] 

NENPQ2.5GPa 

orthorhombic 

0.00 

2 x 10“ J 

48K 

EU 

[Ni(HF 2 )(3-Clpy) 4 ]BF 4 

monoclinic 

0.88 


4.86K 

P»7l 

PbNbVaOs 

tetragonal 

-0.05 

7.8 x 10 - ’ 

104K 

HO 

SrNi 2 V 2 O s 

tetragonal 

-0.04 

8.7 x It) -3 

100K 

HU 

Y 2 BaNi() 6 

tetragonal 

-0.030 

5 x It)* 4 

240K 

[35] 

DTK 

tetragonal 

4.05 

8.2 x It)- 2 

2.2lv 

m 

[Ni(HF 2 )(pyz) a |SbF 6 

tetragonal 

1.22 

1.1 x It)* 1 

OK 


SraNiPtOa 

rhombohcdral 

8.8 


11K 

gd] 


D < 0 in I he D » ./ limit it Ls the M = N - 2 sector that instead determines the 
saturation field, while for intcrmedbite D a discontinuity in M develops at .1 first 
order transition between tlie fully polarized and canted AT AFM pliases.®! 

Gapped and gapless phases respond differently to an applied magnetic field. For 
the longitudinal magnetic field B: — B the Ising AFM, Haldane, and QPM phases 
are all gapped and therefore stay in the M — (I magnetization sector until the field 
reaches some critical value B c . In the ease of the Ising AFM phase, there is a first 
order phase transition the so-called spin flop transition that occurs, whereby 
one sublattice of spins that wils previously aligned against the field “flops" over 
and becomes partially aligned with the field instead, while at the same time the 
oilier sublattice of spins reduces its relative polarization. These spins retain some 
aiitiferromagnctism by developing a staggered magnetization in the XV plane 
the so-called canted AT AFM phase. In the case of the Haldane and QPM phases, 
there Ls instead a continuous phase transition at B c , with field-induced critical 
points belonging to the BEC universality class. 

In Fig. [t] we plot the magnetization curves for = 0.01 for zero-field ground 
states in the Haldane, XV AFM and QPM phases. It Ls easy to see that the gapless 
AT’ AFM phase immediately responds to the applied field, while the magnetization 
curves for the Haldane and QPM phases show signs of a small gap, with B,/B , = 
0.05. All three curves appear similar in their approach to the fully polarized state. 
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Fin. MuguellBttUcui uuil dilfereuliul suMupUbiUty for maguetic field*. applied uluiig the longitu¬ 
dinal aids. For each value of the single-ion anisotropy !>, the Interchain coupling is set in J' = o.oi 
with Invent* temperature /l = L sufficient to reach the ground state regime for L = 32. 


with the increase in slope* due to u combined effect of low spin number (S = 1, 
where we remind the reader that for the most quantum spin .S’ =1/2 the slope 
diverges at B„ as is known from the Bet lie ansaU) ami reduced dimensionality. 
Those features ore most easily seen in t he corresponding differential susceptibility, 
\ = OM/OB. Another feature to notice Is the pronounced hump in the differentia] 
susceptibility around B/B , = 0.6 when D — 0. followed by a local minimum. This 
feature lias previously been pointed out by Kashumikov el of. in t he magnetization 
curve of the Haldane phase of a ID chainand Ls characteristic of the Hiddane 
phase. 

4.5. Low Lying Excitations 

One of the most characteristic features of the Haldane phase is the Haldane gup. 
This arises directly in the Hiddane conjecture due to the difference between integer 
and half-odd-integer spins. It was also one of the first features to be confirmed, with 
early exact diagonalization results showing a finite gap for the spin-1 Heisenberg 
moddPThc accurate determination of this gap was also one of the first applications 

















January 5, 2015 1:30 WSPC/INSTRUCTION FILE 


lmldane-reviow-vS 


14 K. Wierachem ft P. Senguiitu 



7=0.010 
7=0.011 
7=0.012 
7=0.013 
7'=0.014 
7=0.015 
7=0.016 


k.K.n) (O.O.n) (O.x.x) (n.x.re) (x.O.x) (0.0.x) 


Fig. 8. Dispersion us interchain coupling is increased from Huldune phase towards tin* Ndel phase 
ul D = 0. Dashed lino represents the Huldune gup of uncoupled chums (./' = (1). Dut.i shown fbi 
length /. - 32 und inverse temperature d = 2L. Adapted from K. Wierschem and P. Scngupta, 
Phys. Rev. Lett. 112. 247203 (2014). 


of the powerful DMR.CJ mot hod.® 

The low lying states can also he estimated with an upper hound estimator 
25ic/\it > ~i, t hat can ho derived from well known stun rules E3; 1 When the .single- 
mode approximation holds reasonably well, this estimator ciui he quite accurate. 
The spectrum of the Haldane phase in ID Is known from DMRG studiesE£2 and 
in particular .it /.* = rr the spectrum Ls sharply peaked at a single-maguon mode 
(corresponding to the Haldaue gap), with a further gap to the higher multi-inoguou 
modes. Thus, the estimator should he very accurate for k = n. We use this to 
calculate the low lying spectrum of coupled chains with A- = rr along the chain 
direction. Our results are shown in Fig. [Sj where the tendency for the gap to close 
as the interchain coupling Is increased is quite clear, as Is the trend for the inodes 
about k = (jr,jr,»r) to become linear, indicative of the fact that the system is 
approaching the gapless N6el AFM phase. 
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LG. Hidden String Older 


A nonlocal string order captures the hidden symmetry breaking in the Haldane 
phase of the spin-1 HAFM in ID. For two sites i and j in a chain, the string 
correlation function is given by 


c so («' ~ ]) = ~ 


S; exp 


j-i 




in 


k=,+i 



' 0 


In the limit that the distance between sites i and j goes to infinity. Cso becomes 
exactly -4/0 in the AKLT slated while in the Haldane state it 1ms been estimated 
to be 0.3743(1)P 

Although it is not possible to generalize the string order parameter to higher 
dimensional correlations, it is possible to calculate the string correlation function 
in systems of coupled spin chains for spins residing within the same chain. In Fig.[H] 
we plot Cso(L/-i- which as L -+ x should scale exponentially to zero in the 
absence of string order, or to a nonzero value when string order is present. Thus, 
C so (L/2) acts as a string order parameter. For comparison, we also show the scaled 
spin stiffness p % L 2 , which shows that below D, = 0.149(1) the system is in the 
gapped Haldane phase, wliile above this value p. Is finite and we are in the gapless 
X)' AFM phase. It is interesting to note that the string order parameter shows a 
qualitative change in behavior at D,. Below’ D c , Cso(L/2) shows very weak decay 
with system size, and may even saturate to a finite value as i —> oo. By contrast, 
in the XV AFM phase, Cso(i/2) exhibits rapid decay with system size, and will 
scale to zero as L —* oc. While this would seem to indicate that there Is long-range 
string order hi the Haldane phase in QlD geometries, a note of caution must be 
sounded. This is because the string correlation function has been argued to always 
scale exponentially to zero with distance for any nonzero value of the interchain 
coupUngEOHJ Thus, it remains a possibility that the correlation length of string 
order grows rapidly at D,. yet remains finite in the QlD Haldane phase. 


4.7. Symmetry protected topological order 

Due to the shortcomings of the string order parameter in more than one spatial 
dimension, it would be nice to find alternative ways to characterize the Haldane 
phase in QlD geometries. Since the Haldane state in ID is a nontrivial SPT state, 
wo might expect it to remain such in QlD systems. Still, direct evidence fur this 
has been lacking. This may have to do with the difficulty of characterizing SPT 
states in general. One common method Is to look at the low lying entanglement 
spectrum,®^ where an overall degeneracy signals a nontrivial SPT stated" While the 
entanglement spectrum can be easily accessed using powerful DMR.G techniques in 
ID. systems in higher dimension pose a tougher challenge. Additionally, there is 
the possibility that the “cut." used to form a hi part it ion of the system during the 
construction of the reduced density matrix may break an off-site symmetry, and so 
misdiagnose SPT phases protected bv such symmetries. 
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Fin. 9- Fiuilt-siw behavior uf (a) the striug order pwuiuclor C sa (L/2) aud (b) tin* stok'd spin 
stiffness t>,L- uemssthe Haldane to XY AFM phase boundary with aspect ratio R = 6. interchain 
coupling V = 0.01. and Inverse temperature J = L. Adapted from K. Wierscbem and P. Scnguptn, 
Phya. Rev. Lott. 112. 247203 (2014). 


Fortunately, in ID and 2D there exists an alternative way to distinguish be¬ 
tween trivia] and nontrivial SPT states. This is the so-called strange correlator of 
You et al ..® which uses a mapping between the space-time correlations at the spa¬ 
tial boundary of a nontrivial SPT state with a trivial SPT state (or vacuum) and 
(through a Laplace transform) the spatial correlations .it the temporal boundary 
between the time-evolved quantum field theories of the respective trivial .uni non¬ 
trivial SPT states. Since the space-time correlations at the aforementioned spatial 
boundary are known to possess long-range (or. possibly quosi-long-rauge) order in 
one and two dimensions.® a long-range or quasi-long-range strange correlator is 
then a direct measure of SPT order in ID and 2D. while short-range behavior is a 
sign of a trivial product state. In 3D. the possibility of a topologically ordered edge 
state allows for the construction of nontrivial SPT phases with short-range strange 
correlators: however, long-range or quasi-loug-range behavior remains indicat ive of 
SPT order. 

The strange correlator Ls calculated as n "mixed" correlation function between 
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the nontrivial SPT state \ip) and I lie trivial SPT state ji2), 

(«l s;s;. + s>s?. |*> 

W) 

In practice, 11) is taken as a sv linnet tic product state, which Is by definition a trivial 
SPT state. The state |V') can be our putative SPT state in question, and thus a 
(quasi-) long-range strange correlator becomes a direct probe of the nontrivial SPT 
character of the state 1 1 ’ 1 ). It should also be mentioned that this only applies to 
states that are gapped and symmetric, i.e. they are either trivial or nontrivial SPT 
states to begin with. For the case considered here, there is a host of evidence that 
the Q1D Haldane phase is just such a state. The strange correlator has recently 
been implemented using QMC methods.^ 

In Fig. |Ifl| we show the finite-size scaling of a finite-size order parameter 'Vl = 
r, Csc( r ) constructed from the strange correlator. For the QlD Haldane phase, 
we see 1 hat remains finite, signaling a nontrivial SPT state. This is true for buth 

2 D and 3D lattices. By contrast, in the QPM phase we see evidence for exponential 
decay to zero, as expected for a nontrivial SPT state. Thus, we have provided direct 
evidence that the QlD Haldane phase Ls indeed a nontrivial SPT state. 



Csc(r-r') 


5. Discussion 

5.1. Implications of SPT order 

The identification of the quasi-ouc-dimensional Haldane phase as a nontrivial SPT 
state in its own right naturally leads us to consider what sort of topological edge 
modes might be present in this system. Although these edge states are yet to be 
direct ly determined using numerical methods, there exists a solid theoretical frame¬ 
work from which we cmi speculate on their nature. First, consider that the Haldane 
phase in ID supports degenerate spin-1/2 edge states. Next, if we couple Haldane 
chains into an .Y-leg ladder, the edge states will form an overall singlet for N even, 
but retain a degeneracy for jY odd due to Kramers theorem. This explains why 
even leg ladders form trivial SPT states, while odd leg ladders form non-triviid SPT 
states. What happens to the edge states as N -> so? For infinitesimal interchain 
couplings, they must be none other than the spin-1/2 Heisenberg ontiferromaguctic 
chain, with a gapless linear dispersion near the antiferromaguetic ordering wave 
vector. Similarly, the edge states of a three-dimensional QlD Haldane phase corre¬ 
spond to the ground state of the spin-1/2 Heisenberg square lattice antiferromaguct, 
which is likewise gapless. It is natural to assume the edge state spectrum will not 
significantly alter as ,/' Is increased while remaining in the qunsi-oue-dimcusiounl 
Haldane phase. Further, as long as the protecting symmetry of space inversion 
about chain bonds remains intact, these edge states are robust to any symmetry 
conserving perturbations. The quest ion then becomes: How can we measure these 
directly, both in simulation as well as in experiment? In this context, we mention 
that a method has recently been proposed for the experimental probe of fractional 
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Fig. II). Flnite-sixe behavior of the strange order parameter In the Haldane phase (D = II, black 
circles) and the ipiantum paramagnetic phase (D = 2 . red aquaria). Results in the ground state 
limit are obtained by scaling the operator string length an M = 8L.V. Solid lines and filled nymlxils 
are results in :iD witli J' = 0.01. while dashed lines and symbols are in 2D with = 0.02. Adapted 
from K. Wierschi-m uml P. Sengupta. Phys. Rev. Lett. 112. 247203 (2014). 


edge sliil.es iii .S’ = 1 Heisenberg chains.^ 

String order has been predicted to decay exponentially, with correlation length 
C ~ for any finite coupling between eliaiiLsJZUZil If this Ls the case. we 

might expect essentially ID behavior when the length of individual chains L £. 
Since spin-1/2 edge states have been observed in doped Haldane chain compounds 
usiug nun-magnetic dopant ions (which break the spin chains into finite segments), 
wo may surmise that the finite segments Lj opr in such systems obey the relation 
{^dopr) •£ where we use the expectation value to reflect, the distribution of L,<„ (lr 
values due to random placing of the dopant ions (for a doping concentration x. we 
expect {Ldope) = 1/x). In a system such as YjBaNiOr, where the interchain coupling 
is very small (./’ ^ 5 x 10 _4 )P5 the resulting correlation length Ls so large (£ = 
•1 x 10 B ) that almost any finite doping x is sufficient to guarantee that { L,i , tllr ) (. 

Indeed, spin-1/2 edge states are observed in the system Y2BaNii_ x Mg x 05 with 
x = 0.05 and x = (1.1 uF 7 " However, it is interesting to speculate what might occur 
iu Haldane materials with relatively strong ./', such as PbNhVaOs and SrNteVaOs: 
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might there exists a critical doping , below which (Lj ope ) > £? If so. axe spin-1/2 
edge states still detectable? 

5.2. Measurement-based quantum computing 

One exciting possible use of the symmetry protected edge states in Haldane chains 
is for measurement-based quantum computation (MBQO)P^In this quantum com¬ 
puting scheme, an entangled state is initially prepared, after which single qubit 
measurements are performed. These measurements destroy the resource state, so 
this method is also referred to as one-way quantum computing. The AKLT state, 
mid Haldane ground states in general, have been shown to be a possible resource 
state for MBQCE5E3 However, in order to achieve universal quantiun computa¬ 
tion, a two dimensional resource state Is required. In this fashion, the spin-3/2 
AKLT state on the honeycomb lattice has been proposed as a universal resource 
for MBQCigW 1 It is interesting to speculate that any VBS state with protected 
edge modes in two dimensions might prove to be a universal resource state for 
MBQC, in which case the Haldane phase of spin-1 Heisenberg chains weakly cou¬ 
pled into a two dimensional array (as in Fig. [ 2 J 1 might be a way to achieve such 
a state. This would be very useful since AKLT states of higher spins require ad¬ 
ditional Hamiltonian terms beyond the bilinear and biquadratic ones needed for 
spin-1. Spin-1 chains also have the distinction that the AKLT state is adlabatically 
connected to the Haldane phase at the Heisenberg point (i.e. as the biquadratic 
term is tuned to zero, the VBS state remains intact). In higher dimensions, the 
Heisenberg ground state tends towards Neel order, which Is the case on honeycomb 
and square lattices. 

6 . Conclusion 

The remarkable properties of the Haldane phase in spin-1 Heisenberg antiferronmg- 
nets have spurred continued interest in Haldane gap materials and their low-energy 
effective spin models. Vet surprisingly, no accurate determination of the phase di¬ 
agram for these models in experimentally relevant geometries has been completed 
until recently. Here, we have presented some of our recent work on these systems, 
inc luding an accurate phase diagram for .1 square superlattice of weakly interacting 
Heisenberg iUitiferromagnctic spin chains. Further, we have shown the magneti¬ 
zation iunl differential susceptibility curves for longitudinal magnetic fields, and 
calculated the low-lying dispersion along the A‘ ; = 71 plane in reciprocal space. 
We haw also presented direct evidence for the symmetry protected topological 
character of the Haldane phase in quasi-one-dimensional geometries. This includes 
extended string correlations along the weakly coupled chains, as well as a direct 
chiuacterization through the long-range behavior of the strange correlator. These 
calculations lead to a picture of the quasi-one-dimensional Haldane phase as a sym¬ 
metry protected topological phase with an effective gapless spin-1/2 Heisenberg 
iUitiferromagnctic state at its surface. 
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In this brief review, we have endeavored to introduce the Haldane phase within 
the current classification of short-range entangled symmetric states of matter 
namely, the symmetry protected topological states. This exciting development is 
stimulating renewed interest in this field. At the same time, recent developments 
in the growth and chemical synthesis of molecule-based magnets allow for greater 
control over the microscopic properties of magnetic materials. The convergence 
of these two lines of inquiry will he a promising direction for future research. In 
particular, the experimental identification of edge states at the surface of a quantum 
antiferromagnet would potentially realize the first known example of a bosonic 
SPT state. Such a state could have potential uses in measurement-based quant um 
computing. We look forward to many exciting developments in this field in the 
years to come! 
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